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Abstract 

We present a unified framework for designing deterministic monotone polynomial time approxima- 
tion schemes (PTAS's) for a wide class of scheduling problems on uniformly related machines. This 
class includes (among others) minimizing the makespan, maximizing the minimum load, and minimiz- 
ing the £p norm of the machine loads vector. Previously, this kind of result was only known for the 
makespan objective. Monotone algorithms have the property that an increase in the speed of a machine 
cannot decrease the amount of work assigned to it. The key idea of our novel method is to show that for 
goal functions that are sufficiently well-behaved functions of the machine loads, it is possible to compute 
in polynomial time a highly structured nearly optimal schedule. An interesting aspect of our approach 
is that, in contrast to all known approximation schemes, we avoid rounding any job sizes or speeds 
throughout. We can therefore find the exact best structured schedule using a dynamic programming. The 
state space encodes a sufficient amount of information such that no postprocessing is needed, allowing 
an elegant and relatively simple analysis. The monotonicity is a consequence of the fact that we find the 
best schedule in a specific collection of schedules. 

Monotone approximation schemes have an important role in the emerging area of algorithmic mech- 
anism design. In the game-theoretical setting of these scheduling problems there is a social goal, which 
is one of the objective functions that we study. Each machine is controlled by a selfish single-parameter 
agent, where its private information is its cost of processing a unit sized job, which is also the inverse 
of the speed of its machine. Each agent wishes to maximize its own profit, defined as the payment it 
receives from the mechanism minus its cost for processing all jobs assigned to it, and places a bid which 
corresponds to its private information. For each one of the problems, we show that we can calculate pay- 
ments that guarantee truthfulness in an efficient manner Thus, there exists a dominant strategy where 
agents report their true speeds, and we show the existence of a truthful mechanism which can be im- 
plemented in polynomial time, where the social goal is approximated within a factor of 1 + e for every 
£ > 0. 
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1 Introduction 



A major question in algorithmic game theory is how the presence of selfish agents affects the approximability 
of various classic optimization problems lISTl . Specifically, the following research agenda was suggested: 
"to what extent is incentive compatible efficient computation fundamentally less powerful than "classic " 
efficient computation?" (as formulated in ifTSl ). Of particular interest are scheduling problems, where jobs 
are assigned for processing to agents, each controlling one machine, and who have some private information 
regarding their machines |[3Tl l5l [30l[T4]| . In this paper, we consider the case of single-parameter agents with 
scheduling problems on uniformly related machines, which was among the first problems considered in the 
area of algorithmic mechanism design ||5l. The private information of an agent is the cost of processing 
one unit of work, which is also the inverse of the speed of the machine. We provide a negative answer to 
the question raised in [31] for scheduling problems on uniformly related machines, by designing (1 + e)- 
approximation mechanisms for these problems. 

Non-preemptive scheduling problems on m uniformly related machines are defined as follows. We let 
the set of machines be denoted by M = {1, 2, . . . , m}. We are given a set of jobs J = {1, 2, ... , n}, where 
each job j has a positive size pj. The jobs need to be partitioned into m subsets Si, ... , Sm, with Si being 
the subset of jobs assigned to machine i. We let si denote the (actual) speed of machine i, meaning that the 
processing of job j takes ^ time units if j is assigned to machine i. For such a solution (also known as a 
schedule), we let Li = (X^^g^. Pj)/ be the completion time or load of machine i. The work of machine i 
is Wi = Ylij^s Pj ~ ■ '■b^'- '■b^ '■^'■^^ ^^^^ the jobs which are assigned to i. We consider objective 
functions which are functions of the machine loads, Li, L2, . . . , Lm. 

We consider a variety of objective functions (social goals). A well-known objective function is the 
makespan, which is the maximum load. The optimization problem of finding a schedule which minimizes 
the makespan is a basic one Il24l |23l |25l |26l ITSll . The problem of finding a schedule which maximizes 
the minimum load, also known as the cover, is the famous Santa Claus problem on uniformly related ma- 
chines (see e.g. 11221 [32l |2j [8l |20l [HJ |2Ti|). Both these problems are concerned with the optimization of 
the extremum values of the set of machine loads. We will also consider the optimization problem of min- 
imizing X^I^Li fi^i) where / is a well-behaved function. We say that a function / is well-behaved if / 
is a non-negative convex (strictly) monotonically increasing function satisfying the additional property that 
if X < (1 + e)y then f{x) < (1 + 0{l)e)f{y). With regard to the problem of minimizing Yl^i fi^i)^ 
we assume that there is an oracle such that given a rational number x it computes f{x) exactly in constant 
timeQ. The most important example of such a function is f{x) = for p > 1 in which case the problem 
is equivalent to minimizing the £p norm of the vector of machines loads. The optimization goal function of 
minimizing the £2 norm (and the goal of minimizing the £p norm for p > 1) of the vector of completion 
times of the machines has been widely studied (see e.g. iflTl [T3l fJ]). The original motivation was mini- 
mization of the average latency in storage allocation applications (rather than worst-case latency), and the 
problem has additional applications in algorithmic game theory [12|. Bansal and Pruhs ifTOl recently stated: 
"The standard way to compromise between optimizing for the average and optimizing for the worst case is 
to optimize the £p norm, generally for something like p = 2or p = 3." 

The setup of mechanism design for single-parameter agents operating uniformly related machines is as 
follows. Agents present bids to a mechanism, where the bid bi of an agent i is the claimed cost per unit of 

'We can loosen this condition by replacing / with a piecewise-linear continuous convex approximation of / (i.e., the approxi- 
mation is well-behaved as well) without affecting the results. We will assume that / can be computed exactly for simplicity. 
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work of its machine (the inverse of its claimed speed). Based on these bids, the mechanism allocates the 
jobs to the machines and also assigns payments to the agents. We assume that each agent is only interested 
in maximizing its own profit, which is its payment minus its (actual) cost of processing the jobs allocated 
to it. A mechanism is called truthful if reporting their true costs per unit of work is a dominant strategy 
for the agents. That is, this strategy maximizes the profit for each agent, regardless of the strategies of the 
other agents. In the case of single-parameter agents, a well-known necessary and sufficient condition for 
truthfulness is that the allocation algorithm is monotone f5]Sl> that is, the allocation algorithm must have 
the property that if an agent i increases its claimed speed (i.e., decreases its bid) while all other bids are 
unchanged, the work allocated to i does not decrease. More precisely, in such a case there exist simple 
payment functions that can be coupled with the (monotone) allocation algorithm to give a truthful mecha- 
nism. If the allocation algorithm runs in polynomial time, and the payments can be computed in polynomial 
time as well, then the resulting truthful mechanism can be implemented in polynomial time. Thus, for 
single-parameter agents, since the problems are typically strongly NP-hard, the primary goal is to design a 
monotone (polynomial time) approximation algorithm with the smallest possible approximation ratio, and 
to show how the corresponding payments can be computed in polynomial time for its outputs. 

An 7^-approximation algorithm for a minimization problem is a polynomial time algorithm which al- 
ways finds a feasible solution of cost at most TZ times the cost of an optimal solution. An 7^-approximation 
algorithm for a maximization problem is a polynomial time algorithm which always finds a feasible solu- 
tion of value at least ^ times the value of an optimal solution (we use the convention of approximation 
ratios greater than 1 for maximization problems). The infimum value of TZ for which an algorithm is an 
7^-approximation is called the approximation ratio or the performance guarantee of the algorithm. A poly- 
nomial time approximation scheme (PTAS) is a family of approximation algorithms such that the family has 
a (1 + e)-approximation algorithm for any e > (the running time must be polynomial in the input size). 
If the running time is polynomial in ^ as well then the PTAS is in fact an FPTAS (fully polynomial time 
approximation scheme). On the other hand, if the running time is quasi-polynomial (logarithmic factors 
of the input size may appear in the exponent), then the approximation scheme (which is not a PTAS) is a 
quasi-polynomial time approximation scheme (QPTAS). Being strongly NP-hard, the scheduling problems 
studied here cannot have an FPTAS unless P=NP. 

A classic PTAS for these problems generally works by restricting the set of allowable schedules and 
approximating over this set, where the details depend on the specific algorithm and the objective function 
considered. Typically, a chief method of restricting allowed schedules is to do grouping and rounding of 
jobs, where given subsets of jobs are seen as identical, and to treat jobs which are very small compared to 
the work that a machine should receive as arbitrarily divisible (or sand). A number of difficulties arise when 
trying to modify such schemes to satisfy the monotonicity requirement (some of which were partially dealt 
with in the past, see below). It is no longer possible to treat similar jobs as "identical", and their exact sizes 
must be considered. Jobs which are small for the machine which receives them are much more difficult; 
such jobs usually do not affect the approximation ratio but which nevertheless need to be assigned very 
carefully in order to satisfy the monotonicity requirement, since even a very small reduction in the work 
when the machine increases its speed is not allowed. Moreover, it is not known in advance which job is 
small on which machine. 

Dhangwatnotai et al. iTTSl used randomization to construct a monotone PTAS for the three main objec- 
tive functions listed above (makespan, cover, and Ip norm), which combined with an appropriate payment 
function they give, implies a mechanism which is truthful in expectation. That is, given a choice of e > 0, 
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their algorithm for this value of e has an approximation ratio of 1 + e for any realization, but the mono- 
tonicity is proved for the expected works of machines. In this weaker notion of truthfulness, the agents are 
not interested in their actual profits but only in the expected ones, that is, the agents are risk-neutral. For 
example, if an agent earns a profit of M with probability -jj then it sees it as a profit of 1, while a human 
agent would very much be interested in the value of M, and if it is large, it would see it as earning nothing 
at all (rather than earning 1 in expectation). Their approach of dealing with the difficulties above is that 
when a machine receives a job of a given rounded size, the actual job is chosen uniformly at random from 
the set of jobs of this rounded size, so the "sizes" of jobs (the expected sizes) are easier to deal with. For 
jobs that are small, a fractional assignment is found (and rounded using randomization). They also derived 
deterministic monotone QPTAS's for minimizing the maximum load and the Ip norm of the loads. A fully 
deterministic (and hence universally truthful) monotone PTAS for minimizing the makespan was given by 
Christodoulou and Kovacs 115 ]. They assign jobs that have almost the same size (are in the same group) 
very carefully in a fixed order (sorted by size) to the machines (where machines are given in a fixed order of 
their speeds). Moreover, they begin by rounding speeds to powers of 1 + e, and round the job sizes to powers 
of 1 + 5 for some 5 e. This ensures that when a speed changes, this change is always relatively large 
compared to the job classification, so the rounding errors introduced by small jobs are not large compared 
to the required change in the work. The authors give a long and technical proof to show that it is possible to 
combine these main ideas and give a deterministic monotone assignment. This approach can be used only 
for minimizing the makespan, since in the scheme of [ 15], machines of similar speeds should either receive 
almost the same work (implied by the makespan), or no small jobs at all, unless no small jobs remain. 
Informally, the small jobs are pushed to the fastest machines. This approach does not seem to work even 
for the similar problem of maximizing the cover, but applying the methods of [15 ] leads to a deterministic 
monotone (2 + e) -approximation for this last objective, given by Christodoulou, Kovacs, and van Stee fl6l 
(the problem was also studied in ll2ni ). 

What can be seen from these previous results is that satisfying the monotonicity requirement would 
become easier if we could simply avoid the notion of small jobs. Then we could calculate with exact job 
sizes (and thus exact loads) throughout. An important contribution of this paper is to show that for any 
given schedule, a highly structured schedule exists, where the ratio of job sizes assigned to a machine is 
unbounded but the jobs types assigned to this machine are restricted in the sense that these jobs are grouped 
into a sufficiently small number of classes. This overcomes the difficulty that it does not seem to be possible 
to actually bound the size ratio of jobs assigned to a machine, but still we would like to use dynamic pro- 
gramming without introducing a notion of small jobs or inexact calculations. The set of possible outcomes 
is independent of the possible speeds, which assists in dealing with speed changes, and finally, the work 
of each machine is very close to its work in the given (original) schedule, which keeps the approximation 
ratio close to 1. This allows us to deal with all of the objective functions mentioned above at once using 
a dynamic programming formulation implemented by a layered graph, having one layer for each machine. 
Unlike previous approximation schemes which use such graphs, a path in the graph corresponds to one spe- 
cific schedule (not to a class of schedules, or a schedule for a set of rounded jobs), and the cost of the path 
(with respect to a goal function) is precisely the cost of the corresponding schedule and not its approximated 
value. That is, there is no rounding or imprecise calculation with respect to relatively small jobs (or any 
other jobs). This makes proving monotonicity much more straightforward, and even simplifies the proof 
of the approximation ratio, and the presentation of the algorithm, compared to previous (non-monotone) 
PTAS's. Our construction works in the same way for all inputs and all objectives, and does not require any 
special cases. Hence we streamline the monotone PTAS for minimizing the makespan [[15 J . Moreover, we 
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provide the first deterministic monotone PTAS's for maximizing the minimum load and minimizing the ip 
norm, which are our main contributions. 

Other related work. For a fixed (constant) number of machines, scheduling problems typically have an 
FPTAS |27j;9nl9||, and even a (deterministic) monotone one for makespan minimization and for maximizing 
the minimum load ll3j|211. The QPTAS of |[T8l for minimizing the Ip norm is in particular a PTAS for fixed 
values of m. Prior to the monotone FPTAS of Andelman, Azar, and Sorani 131 for makespan minimization, 
Auletta et al. 161 gave the first deterministic monotone algorithm for this problem (where the number of 
machines is fixed), with an approximation ratio of 4 + e. 

In what follows we discuss the case where the number of machines is part of the input. It was shown 
by Hochbaum and Shmoys that the makespan minimization problem has a PTAS for identical (equal speed) 
machines |[25l and for uniformly related machines ||26l . All optimization problems studied here, including 
maximizing the minimum load and minimizing the £p norm, are known to have a PTAS for identical ma- 
chines ll32l [Tll2]|. and for uniformly related machines 1111201. As for monotone algorithms for the makespan 
minimization problem, before the papers ifTSlfTSl mentioned above. Archer and Tardos [5| gave a random- 
ized 3-approximation mechanism for minimizing the makespan which is truthful in expectation only. The 
ratio was later improved to 2 H| (and eventually to 1 + e |[T8l ). A deterministic monotone algorithm of 
approximation ratio at most 5 was given in Q, and Kovacs improved the ratio to 3 and then to 2.8 Il28ll29l . 

Proof overview. Our proof consists of two parts. In the first one, we define several properties which a 
structured schedule should have, and show that every schedule has a similar schedule which has such prop- 
erties. As stated above, similarity is measured by allowing only a very small change in the work of every 
machine. For the proof we introduce a notion of a fractional schedule, where some (relatively small) jobs 
may be spUt over multiple machines. For any (integral or fractional) schedule, we can define a magnitude 
vector with a component for every machine. Unlike previous work, where the magnitude of a machine corre- 
sponded directly to its work (or the largest job assigned to it), we use the magnitude component of a machine 
as an upper bound for the size of any job which is assigned to it, but if a component of the magnitude vector 
is different from the previous one, we require that the value of this component matches (approximately) 
the work of the corresponding machine. There are several ways to define a magnitude vector for a given 
schedule. A possible solution to the dynamic programming can be viewed as a process where we create 
the magnitude vector component by component (for a list of machines sorted by non-decreasing speed); 
increase the magnitude of the current machine (as opposed to keeping the same magnitude of the previous 
machine) only if keeping the same magnitude as for the previous machine would result in a violation of the 
upper bound on the maximum size of any job assigned to the current machine. This novel approach allows 
additional flexibility in the set of allowed schedules. 

For a given integral schedule, where the works of the machines are increasing with the speeds, we show 
that a fractional schedule exists where the total size of very small jobs which are (partially) assigned to 
machines with high work is small, and the work on each machine is the same as the work in the integral 
schedule. We then refine this result by constructing an integral schedule where no very small jobs are 
assigned to machines with high work, the works of the machines are all close to the original works, and an 
additional technical property holds. However, despite the works being close to the original works, they may 
no longer be sorted in the resulting schedule (though if the works of two consecutive machines are unsorted, 
then the difference between their works is very small). Searching for unsorted schedules causes technical 
difficulties for the algorithm which should find a structured schedule, while a postprocessing step of sorting 
may harm monotonicity. We therefore do one extra step to create a final integral schedule in which the works 
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are sorted again (but still very close to the original works) and several structural properties hold. We do not 
use rounding, but jobs are partitioned into mega-classes and mini-classes according to their size, and we 
apply re-assignment of jobs in every class to comply with the required structure. For a given schedule, some 
classes of jobs can turn out to be too large for some machines, while they are very small compared to the 
work of other machines. These jobs are combined into chunks called "alternative jobs". Since this process 
can be applied in particular for an optimal schedule (for each one of the studied problems), there exists a 
schedule where works are very close to the works in an optimal schedule, and the structured schedule has 
an objective value which is close to optimum. 

Once we show the existence of such a schedule, we can turn to the design of an algorithm which finds it. 
We use a dynamic programming formulation which is based on the structural properties. By the structural 
properties and the existence of a magnitude vector, it is only necessary to have a small number of components 
of this vector in the state space. A preprocessing step is performed, where all possible types of alternative 
jobs are created. While a job will belong to a number of sets of alternative jobs, every solution will use it at 
most once as a part of an alternative job (or possibly it will simply be assigned as a job). Thus, we find an 
optimal solution out of a given class using a polynomial time algorithm, and this optimal schedule is then 
guaranteed to be close to an overall optimal schedule, as well as being monotone. 

2 Preliminaries 

For our results, we let e be a small constant such that < e < ^ and ^ is an integer power of 2 denoted by 
r > 5 (i.e., e = ^). Throughout the paper, for a solution A we denote by A both the solution and the value 
of the objective function for this solution. Without loss of generality, we assume that < pi < p2 < ■ ■ ■ < 

Pn- 

An integral schedule is a function S : J ^ M. We let Wf = J2jeJ-S{j)=iPj ^^^^^ work of 

machine i in the integral schedule S). A fractional schedule is a function X : J x M — > [0, 1]. The value 
X{j, i) is the fraction of job j assigned to machine i, and the following condition (that every job is assigned 
completely) must be satisfied: 
(Fl) For every j G J, J2ieM = 1- 

Let W/^ = Pj ' ^(ii the total fractional size of jobs of machine i, and let = 2°^ , where 

af = |"log2 W^] , be its rounded value (if = then af = -oo and = 0). We call the work 
of machine i'mX (as for integral schedules) and is the rounded work (also for integral schedules). A 
fractional schedule is valid if it satisfies condition (F2): 

(F2) There is a partition J = Jz{X) U Jr(X) (Jz(X) n Jr(X) = 0), such that if j G Ji{X) then there is 
a unique value i e M such that X{j, i) > (and therefore X{j, i) = 1), and if j G Jr{X) and X{j, i) > 
then Pj < eW^. 

Note that the partition in (F2) is not necessarily uniquely defined. Every integral schedule S induces 
a valid fractional schedule X with the same jobs assigned to every machine as follows: let i) = 1 if 
S'(j) = i, else i) = 0. Furthermore, we let Jr{X) = {j e J : pj < and Jz{X) = J\Jr{X). 

Note that Wf = for i = 1 , . . . , m. X is called the (valid) fractional schedule induced by S. On the 
other hand, every valid fractional schedule X for which X{j,i) G {0, 1} for all j ^ J,i & M induces 
an integral schedule S with the same works by setting S{j) = i for the value of i for which X{j,i) = 1 



5 



(this value of i is unique due to (Fl)). S is called the integral schedule induced by X. In what follows we 
use the term schedule for an integral schedule. We let Lf = be the load of machine i in the schedule 
S. The first part of Claim [T] follows from an observation in (TH], and it is easy to show the second part 
(see Appendix lA.ll ). For all cases, we conclude that if machines are sorted by non-decreasing speed, it is 
sufficient to consider optimal schedules where the works are non-decreasing (as a function of the indices). 

Claim 1 Assume that si < S2 < • • • < Sm- There exists an optimal schedule S for the problem of minimiz- 
ing YllLi f{Li) where f is a well-behaved function, which satisfies Wf < W2 < • • • < W^. There exists 
an optimal schedule Si for the makespan minimization problem which satisfies W^^ < VF^^ < • • • < W^. 
There exists an optimal schedule S2for the machine covering problem which satisfies ^ < W2 ^ < • • • < 

" m • 

3 The existence of near-optimal highly structured solutions 

We define a partition of J into mega-classes. For k £ Z, let I^, = (2^^ , 2*^+^], and let mega-class k be 
{j € J : pj G Zk}- We say that an integer k dominates the integer k' if k > k' + r. Mega-class k dominates 
mega-class k' if k dominates k'. If j, j' belong to mega-classes k, k', respectively, such that mega-class k 
dominates mega-class k', then pji < epj. This holds because pj > 2^ > 2^ +r+i _ i . 2*^'+^ > 1 • p-,^ 
since k' + 1 < k — r and e = 2^"^. We refine this partition and consider the partition of J into mini- 
classes as follows. Denote hy K C Z the set of indices of non-empty mega-classes (clearly \K\ < n). Let 
A = [logi+^ 2]. For A; G K and < £ < A - 1, let = {2^ ■ (1 + e)^ 2^= • (1 + eY+^]. The mini-class 
{k, tj is the set of jobs of mega-class k whose size is in Ik,i. Note that (1 + e) > (1 + e)'°si+^' ^ ^ 2 

and thus the partition of J into the mini-classes is a refined partition of the partition into the mega-classes. 

Given a set of consecutive mega-classes ki, . . . ,k2 where k2 > ki, with the job set J consisting of all 
jobs of J with size in the interval (2'^^ , 2^^2+1]^ and letting g = 2^^^, we create an alternative set of jobs 
that will possibly replace J. These alternative jobs have size in the interval {g, 2g] (except perhaps for one 
alternative job that may be smaller). To create these alternative jobs we partition J into subsets each of 
which has total size at most 2g such that no two subsets can be united keeping this condition. A set of 
subsets satisfying this condition has at most one subset whose total size is at most g. We create these subsets 
by picking in each step a maximal prefix of the jobs in J (where J is sorted according to the indices of 
the jobs, i.e., by non-decreasing size) with total size at most 2g and remove the selected jobs from J. This 
algorithm is equivalent to applying the bin packing algorithm Next-Fit Increasing (NFl) using "bins" of size 
2g; once a subset of total size at most g is picked, all further subsets (if any exist) have total sizes above g. 
The algorithm sometimes decides to replace J with the alternative jobs, and in this case we partition these 
alternative jobs into separate mini-classes which we call alternative mini-classes. The alternative mini-class 
{k, €) contains all the alternative jobs with size in 1^/, resulting in at most A + 1 alternative mini-classes. If 
the algorithm decides to replace J with alternative jobs, then in the output of the algorithm each alternative 
job is replaced with the original jobs which were combined to form it, and this is done just before returning 
the output (the work of each machine is not affected by this change). Since there are at most n non-empty 
mega-classes, there are O(n^) different sets J that possibly the algorithm replaces with alternative jobs. 
Thus creating all the sets of alternative jobs takes 0{iv'). Note that one job can be contained in multiple 
alternative jobs, but at most one of these alternative jobs will be used. 

Definition 2 An integral schedule respects the alternative jobs of mega-classes ki, . . . ,k2, where k2 > ki, 
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if every pair of jobs j, j' with size in the interval (2*=i, 2^=2+1] which are within a common subset ( that is, 
should be combined into one alternative job with possibly other jobs), are scheduled on a common machine. 

The motivation for this definition is that these jobs can be easily replaced by the alternative job to which 
they belong without affecting the works of the machines. 

Definition 3 A vector a = (ao, ai, . . . , am) (of length m + I) whose components belong to {—00} is 
called a magnitude vector if ao = —00, for i = 0,1, ... ,m — 1, ai < Oj+i and if a^ 7^ a^+i then a^+i 
dominates ai (i.e., Oj+i > ai + r + 1). 

We now define the signature vector 6 of a magnitude vector a. The number of components in b is the 
number of distinct values among the components of a excluding oq, denoted by r(a). Each component 
t = 1,2, . . . , r(a) of 6 is a pair bt = (Ct, I't) such that ^1 = 1, and for 1 < t < r(a) and < ^ < ^t+i — 1 
(where ^T-(a)+i = m + 1) we have at = vt. That is, the value is always the first machine which has a 
larger component of a than the previous machine and this component is vt. For every t = 1,2, . . . , T{a) — 1, 
we let J\a) = {j £ J : 2"'+''+'^ < pj < 2'^*+!-''}. 

Observation 4 For every job j and every magnitude vector a with its signature vector b, there are at most 
two values of t £ {1, ... , T{a)} for which pj G {2'^^~^, 2^*+''+^], and if there exists at least one such value 
oft, thenj ^ U0j^{a). 

Proof. By the definitions above, for every j there are at most two values of t for which pj G (2'^*"'', 
2i^t+r+ij (since U0+1 > + r + 1 for every 9). Moreover, if G (2"^-'', 2'"'+'"+i], then for every 6 < t, 
we have j ^ J^ia) because pj > 2'^*~^ > 2'^'>+^~^', and thus j is too large to be in J^{a). \i 9 > t, then 
Pj < 2^*+^+i < 2'^»+''+\ and thus j is too small to be in J^(a). ■ 

Definition 5 A valid fractional schedule X is consistent with a magnitude vector a if 1) for every job j 
and machine i, if X(j,i) > then pj < 2''»+''+^, that is, machine i does not contain parts of jobs of a 
mega-class higher than Ui + r, and 2) if ai ^ ai-i (for i G M) then ai = otf(= [log2 W/^~\ ). 

Observation 6 If a valid fractional schedule X is consistent with a magnitude vector a and < < 
■ ■ ■ < then for every i G M, we have ai < af . 

Definition 7 A pair {X, a), where X is a valid fractional schedule, and a is a magnitude vector such that 
X is consistent with a is called favorable if for t = 1,2, . . . , r(a) — 3, we have 

m 

J2 E Pj-x{j,i)<2^^^^+^+K 

This condition ensures in particular that the total size of parts of jobs whose mega-class is dominated by 
mega-class fj, assigned to a machine of index at least ^^+3, is relatively small compared to the work of that 
machine. This holds since 2'^*+i+''+i < 2''*+2 < 2''*+^-'' = e ■ 2"*+^ = e ■ Wf^ . 

We define several processes in which a vahd fractional schedule is modified into a different vaUd frac- 
tional schedule. These processes are defined algorithmically but they are not a part of the final algorithm, 
but only of the proof that a highly structured integral schedule must exist. 
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FNFI. For a subset of jobs J' C J and a set of bounds f/i, . . . , Um (for the m machines) such that 
Ylij&j'Pj — S2=i '^he Fractional Next-Fit Increasing (FNFl) algorithm creates a fractional allocation of 
these jobs in the following way. Let i = 1 be the first active machine, and for every j G J' let qj = pj. In 
every step, fnfi picks the minimum index job j G J'. It allocates /3 = mm{qj, [/, } processing time of this 
job to machine i. It decreases both Ui and qj by /3. If Ui = 0, then it increases i by 1, and if qj = 0, then it 
removes j from J'. FNFI repeats this step until i = m + 1 (and J' = must hold, these two events happen 
simultaneously since J2jeJ' Pi ~ Si^i ^*)- ^^^^ sometimes used to reassign a subset of jobs in a valid 
fractional schedule, so that the total sizes of jobs of this subset assigned to each machine is unchanged (i.e., 
the bounds Ui are given by assignment of the jobs of the subset in the original valid fractional schedule). 
This is done only in situations where it is ensured that the resulting fractional schedule is valid. 

Definition 8 A valid fractional schedule X is compatible with FNFI ;/ running FNFI on the input job set 
Jr(X) with the set of bounds Ui, . . . , Um such that U = Y2j£jji{x) Pj-^U^ ^) allocates exactly pj ■ X{j, i) 
time units of job j to machine i for every j € J^{X) and all i G M, that is, keeps the valid fractional 
schedule unchanged. 

Round-FNFI. On several occasions, given a valid fractional schedule X, which is compatible with FNFI, 
we will apply the following rounding procedure, called Round-FNFI. Assign each job j € Jr(X) completely 
to the minimum index i such that X{j,i) > 0. Since in the assignment process of FNFI each machine 
receives at most two jobs which are not completely assigned to it, the one of smallest index and the one 
of largest index, the resulting fractional schedule induces an integral schedule S in which each machine 
may have additional parts of at most one job (the one of the largest index assigned to this machine by 
FNFI), and may have less parts of at most one job (the one of the smallest index assigned to this machine 
by FNFI). Since by condition (F2) each fractional job j G Jr{X) on machine i (that is, every j G Jr{X) 
such that X{j,i) > 0) has size pj < eW^ < 2£W^ , we conclude that for every i G M we have 
(1 — 2e)W^ < Wf < (1 + 2£)W^ . We say that the integral schedule S is created by applying Round- 
FNFI on X. 

Lemma 9 Given a schedule S : J ^ M such that Wi < W2 < ■ ■ ■ < there exists a favorable pair 
{X, a) where = for i = 1,2, ... ,m, and X is compatible with FNFI. 

Proof. First, as described in Section |2l S induces a valid fractional schedule, here denoted by Xs, with 
the same sequence of works. For i = 1,2, ... ,m, let qi = af. Since Wf < W2 < ••• < W^, 
we have q^ < (/j+i for alH = 1,2, . . . ,m — 1. We define a magnitude vector = (og , af, . . . , a^) 
as follows. We let Oq = —00, for i = 1,2, ... ,m, if qi < af_^ + r, then af = af_'^, and otherwise 
af = qi. The valid fractional schedule Xs is consistent with since for every i ^ M either af = qi or 
Qi < a'i-i + r = af + r. In both cases, the size of any job assigned (completely) to machine i cannot exceed 

We next consider the nonempty set of pairs {X' , a) such that X' is consistent with a, and such that for 
i = 1,2, ... ,m, = and a^ < (the set is indeed nonempty by the existence of {Xs, a^)). 

Among all the possible choices for X' and a, we consider one such that the vector a has a signature vector 
with the smallest number of components, and (as a secondary objective, i.e., among such solutions which 
minimize the number of components in the signature vector) | J]r(X')| is maximized. Based on X' we will 
define X (by applying FNFI on Jr{X')), and X will be shown to be a valid fractional schedule satisfying 
the lemma. 
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We modify X' by reassigning the jobs of Jr{X') using FNFI with the set of bounds Ui, . . . , Um such 
that Ui = Pj^'U: 0- We denote the resulting fractional schedule which is compatible with FNFI 

by X. We argue that X satisfies (F2). We define Jr{X) = Jk{X') and show that if j G Jr{X) and 
i G M satisfy that X{j,i) > 0, then pj < eW^ . Since the works of the machines are sorted in a non- 
decreasing order, it suffices to show that for j G Jr{X) and i such that X{j, i) > 0, there exists j' > j, 
j' G Jk(X') and i' < i such that X'{j', i') > 0, since in such a case pj < pj/ < eW^ < sW^ . Assume by 
contradiction that this claim does not hold for j and i. Then, since FNFI assigns job j (possibly partially) to 

machine i, Y.j'eMX'):j'<jPr < T!-y=i U-/, however T,j,(=j^(^x'y.j'<j Pj' > E7=i since no other jobs 
of Jk(X') are assigned by X' to the first i machines. Therefore X is indeed a valid fractional schedule. 

We claim that X is consistent with a. It suffices to prove that in every prefix of machines 1,2, ... ,i, the 
maximum size of a job j such that 7) > for some 1 < 7 < i does not increase when we replace 
X' by X. Let j be a job of maximum size which is assigned in X (possibly fractionally) to a machine 
7 G {1,2, If j G Jz{X) = Jz(X')then X'(j,-f) = X{j,-/) = 1, and the claim holds. Otherwise, 

j G Jr(X). There exists / G Jr(X') and < 7 such that X'{j', i') > and / > j as we showed above, 
and the claim holds as well. 

Last, we prove that {X, a) is a favorable pair. Let t be such that 1 < t < T{a) — 3. Let j G J be such 
that there is i G [^4+3, m] with X{j,i) > and pj < 2*^*"^. If there is no such job, then we are done. 
We have j G Jr{X) = Jr{X') because W^^ > 2"^ > 2"*+-^ > 2"' • 2^^ = ■ 2"' > j^pj where the 
first inequality holds by Observation |6j so if j ^ Jr{X) then X{j,i) = 1 and we can add j to Jr(X), 
contradicting our choice of X' . Consider the machines ^f+i = {■^t+i, • • • ,it+2 — 1}- If all jobs assigned 
(possibly fractionally) by X to these machines have size of at most 2'^*+''+^, then we can redefine Oj' for 
i' G AtJ^i to be ft, contradicting the minimality of length of the signature vector of a. Consider a job j' 
such that there is i' G At+i for which X{j',i') > and pj' > 2''*+''+^. By the existence of j G Jr{X) 
with size at most 2''*"^, such that a part of it is allocated to a machine of higher index, we conclude that 
/ G Jz{X) since X is compatible with FNFI. We also have pj' < 2"^+^+''+^ < 2"^+'^-''-^ < e ■ ^. 
If ET=5*+3 Ej":r,,,<2''.-r Pj" ■ X{f,^) > 2-*+i+^+i, then ^7=^^^, E,":p^,„<2n-. Pj" ■ X{j" > p-,. 
In this case, we add j' to Jr{X), and modify X as follows. We consider a replacement of the position 
of j' with the position of a set of fractions of jobs (where each such job has size at most 2''*"'' = e2'^* < 
g22!^t+i~i < £2^^ ^ e^W^ for every 7 G At+i, and belongs to Jr{X')) of total size py which were 
previously assigned to machines with index at least ^4+3. The resulting schedule indeed satisfies (F2) since 
the jobs which take the place of j' are smaller than e^W^ for every 7 G ^t+i while pji < e ■ W^^,^. Thus, 
the resulting valid fractional schedule is consistent with a, contradicting our choice of X' since | Jik(X')| is 
not maximal among valid fractional schedules consistent with a (and having the required properties). ■ 

Definition 10 A schedule S is almost consistent with a magnitude vector a if for every i = 1,2, . . . ,m, the 
set of parts of jobs assigned to machine i does not contain any part of a job of a mega-class higher than 
Oj + r, and if at 7^ aj_i (for i G M) then \ai — af | < 1. 

Definition 11 A schedule S : J ^ M is good if the following properties hold. 

1. There exists a magnitude vector a such that S is almost consistent with a, and furthermore for every 
t = 1,2, . . . , T{a) — 4 there is no j and i > (^t+A such that pj < 2''*^'' and S{j) = i. 

2. For every t=l,2,..., r(a) - 1 if J\a) = {j e J : 2"'+''+'^ < pj < 2''^+^-''} / 0, then X respects 
the alternative jobs of mega-classes vt + r + 1, . . . , vt+i — r — 1. 
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Lemma 12 Given a schedule S : J ^ M such that < W2 < ■ ■ ■ < t^^^^ exists a good schedule 
S' : J ^ M such that for i = 1,2, . . . ,m, we have 

(1 - 12e) • Wf < Wf < (1 + 12e) • Wf . (1) 

Proof. By Lemma |9j there exists a favorable pair {X, a) where W^^ = Wf for i = 1,2,..., m, and X 
is compatible with FNFI. First, for every t = 4, 5, ... , T(a), we reschedule all parts of jobs j such that 
Pj < 2^^-'-^~^ and for which there exists i > such that X{j, i) > by moving them to machine We 
denote by X the resulting fractional schedule. We next bound the value of in terms of for every 
i G M. The work of i may increase (if i = £,t for some t = 4, 5, ... , r(a)). Since {X, a) is a favorable 
pair, the amount of this increase is at most e ■ < 2eW^ , since 2''*-2+r+i ^ 2'^*~^ = eW^ . Next, 
we bound the total size of parts of jobs removed from machine i (for 2 < i < m). Let t' be the maximum 
index such that < i (which must exist since ^1 = 1). Then, for every t = 4, 5, . . . , t', we may have 
removed a total size of at most 2''*-2+''+i < f • 2*"* from machine i (and move these parts of jobs to machine 
it)- Thus - < f • Ya^a 2^"* < e • 2^"*' < 2e • . We conclude that for every i, we have 
(1 - 2e)W^ <W^ <{l + 2e)W^. 

Let Jk(X) = Jk(X). We observe that X is a valid fractional schedule which is compatible with FNFI 
(similarly to the bounds on such jobs in Lemma|9j it can be shown that if a job moved to machine i, then its 
size is below eW^ , since 1 — 2e > e). We now apply Round-FNFI on X to create an integral schedule S. 
Every j G Jir(X) such that X{j, i) > has size pj < eW^ < 2eW^ , so for every i ^ M we have 

(1 - 4e)W/^ <Wf <{l + A£)W^ . (2) 

The maximum size of a job in a prefix of machines in S is the same as in X, and a job moved from its 
position in X to a new position on machine in X has size at most 2^^*^^ < e2^* = e2°-it . 

Consider the set of jobs J*{a). Since X is consistent with a, for every j € and i < ^t+i, we have 

X{j, i) = 0, and since the maximum size of a job in a prefix of machines did not change, S{j) > i. Since 
W^^^^ = 2'"*+!, we have for all j G J*(a) and i > ^t+i thatpj < 2''^+^-'' = e-W^^^^ < e-W^. We remove 
the jobs in J*(a) from their positions in S, and we will schedule the alternative jobs instead (which gives a 
schedule of the original jobs which respects the alternative jobs of mega-classes vt + r + l, . . . , vt+i —r — 1). 
For every i G M we let Ui be the total size of jobs in J*(a) which are assigned to machine i by S. The set 
of machines i for which f/j / is contained in the interval [it+i,it+4\ where if t + 4 > T(a), then we let 
it+i = m. We apply fnfi to fractionally schedule the alternative jobs, followed by Round-FNFI. This is 
done for every value of t for which J* (a) ^ sequentially. We denote by 5' the resulting integral solution. 
Let i G M. There are at most four values of t for which i participated in the process of the rescheduling 
of J^{a). As a result of applying Round-FNFI for the alternative jobs for all t, every machine i can have 
at most four additional parts of jobs and less parts of at most four jobs, all of which have size of at most 
eWi^ < 2eW^. Thus, - ^eW^ < Wf <Wf + 8eW,^. Using ©, we get ©. 

The integral schedule S' is almost consistent with the magnitude vector a. To see this claim, first 
observe that no job is too large: if the maximum size of a job on machine i in S' is not the same as in 
S, this maximum size job j G J*(a) is moved from its position in 5 to a new position on machine i, and 
therefore pj < 2'^*+i~'' = e2'^*+^, and > a^^^-^ = vt+i- The claim holds because for every i, we have 
|af - af I < 1 since 1 + 12e < 2 and < 2. ■ 
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Definition 13 A schedule S is quasi-consistent with a magnitude vector a if for every i = 1,2, ... ,m such 
that it ^ i < ^t+i, the set of jobs assigned to machine i does not contain any job of a mega-class higher 



Definition 14 A schedule S : J ^ M is structured if the following properties hold. 

1. There exists a magnitude vector a such that S is quasi-consistent with a, and furthermore for every 
t = 1,2, ... , r(a) — 5 there is no j and i > ^t^^ such that pj < 2''*"'" and S{j) = i. 

2. For every t = 1,2, . . . , r(a) — 1, if J^{a) ^ 0, then S respects the alternative jobs of mega-classes 
vt+r + 1,.. .,1't+i - r-1. 



4. For each pair of jobs j, j' ^ UtJ*(a) belonging to a common mini-class, if j < j', then S{j) < S{j'). 

5. For each pair of alternative jobs resulting from the set J^{a) belonging to a common alternative 
mini-class such that the size of j is smaller than the size of j', the following holds. If S schedules the 
original jobs in j and j' on machines i and i', respectively, then i < i'. 

Tlieorem 15 Given a schedule S : J ^ M such that Wf < < • • • < there exists a structured 
schedule S* : J ^ M such that for i = 1,2, . . . ,m, we have 



Proof. Let S' be the good schedule that is based on 5 as established in Lemma [12] We apply a sorting 
procedure of the works of the machines similarly to the one of ifTSll . In this procedure we are given as 
an input a partition of the jobs into subsets Ji, . . . , Jm, we create a new partition of the jobs as follows. 
For i = l,2,...,m — Iwe assume that we are given the subsets Ji, . . . Jm and we choose the set of 
jobs scheduled on machine i (possibly modifying the remaining subsets). For each mini-class (including 
the alternative mini-classes), we temporarily replace the jobs m Jy (i' = i,i + 1, . . . ,m) from this mini- 
class with the smallest set of jobs of this mini-class which are still available (i.e., they are not scheduled 
on machines with indices smaller than i). We pick the set of jobs which have minimum total size as the 
set Ji, possibly swapping locations of jobs in the same mini-class. Note that due to our use of alternative 
mini-classes, the jobs that are inside these alternative jobs might not be allocated in order of their size (but 
still in a fixed order according to the size of the alternative jobs). Consider a pair of consecutive machines 
i,i + l, then the resulting work of machine i is not larger than the resulting work of machine i + 1, since the 
set of jobs allocated to machine i + 1 were available for allocation to this subset of jobs when we picked J 
for machine i (the jobs taken by J are replaced by other jobs of the same mini-class when we choose J+i, 
which cannot be smaller). 

We apply the sorting procedure on the partition defined by S'. The output of this procedure is an integral 
schedule denoted by S*. Clearly, Wf < W2'' < ••• < W^* . Moreover, properties |2l S] and [5] in the 
definition of structured schedules are satisfied. We next prove Q for every machine i. Every machine 
i receives a subset of jobs which is based on a subset of jobs allocated to some machine i' in S', after 
swapping pairs of jobs within a common mini-class. Therefore, 



than vtJ^i + r, and if ai 7^ aj_i (for i G M) then \ai — a'. 




3. Wf <wi <■■■ < W^. 



(1 - 14e) -Wt <Wt < (1 + 14e) • 



(3) 




< Wf <{l + e)-W^ 



(4) 



1 + e 
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Fix a machine index i. When we choose the set J'i for machine i, at least one of the original i subsets 
J'l, . . . Ji (up to swapping some locations of pairs of jobs within common mini-classes and alternative 
mini-classes) remains available. Let i" denote its index. The total size of the jobs in this subset is at most 
W^^' • (1 + e) < VF^» • (1 + 12e)(l + e) < Wf ■ (1 + 12e)(l + e) where the first inequality holds by Lemma 
[T2I and the second inequality holds by the monotonicity of works in S. Therefore, machine i receives in S* 
a total work of at most Wf ■ (1 + 14e). We next prove the other inequality, that is, (1 — 14e) • Wf < Wf. 
The schedule S has at most i — 1 machines which receive work strictly below Wf. Therefore, in S' there 
are at most i — 1 machines which receive work strictly below Wf ■ (1 — 12e). In S* the number of machines 

with work strictly smaller than ^' fi'^"^^^'' cannot exceed i — and due to the monotonicity of the works in 
S*, the claim holds. 

Finally, we prove property [J of structured schedules. The integral schedule S* is quasi-consistent with 
the magnitude vector a. To see this claim, consider a machine i such that < i < ^t+i and denote by j 
the maximum sized job on machine i according to S*. We need to prove that pj < 2'^*+i+''+i. The set of 
jobs J'i which the sorting procedure allocated to machine i was scheduled on a machine i' in S' (possibly 
swapping pairs of jobs in common mini-classes). A job j' of the same mini-class as j was allocated to 

machine i' in S'. Recall that < Wf < (1 + e) • W^' . By Lemma [T2l S' is almost consistent with a, 
and therefore pji < 2"*'+^+^. Since j and / belong to a common mini-class, they also belong to a common 
mega-class, and thus we also have pj < 2"*'+''+^. In order to prove that pj < 2^^*+!+^+^ it suffices to show 
that i' < Assume by contradiction that i' > We have > W^^ ^ > 2"*+'^^^ which holds since 

the works in S are monotonically non-decreasing and 2W^^^^ > W^^^^ = 2 ^'+2 = 2*^*+^. On the other 

hand, Wf < W^^^^ < 2"'+^ By (© and Lemma[l2l we have Wf > ^ > '^'^l+^Z^^'' > W^^{1 - Ue). 
Therefore, we get 



> 



Wf{l-Ue) 

1 + 14e 

1 - 14e 



> 2''*+2-i 



1 + 14e 



> 2^*+2-3 > 2'^*+i+'-2 , since e< — 

32 

contradicting r > 5. Therefore, S* is quasi-consistent with the magnitude vector a because for every i, we 
have |af — af * | < 1 since 1 + 14e < 2 and j^Yii < 

Fix a value of t = 1,2,..., r(a) — 5, it remains to prove that there is no j and i > such that 
Pj < 2*^*"^' and S* schedules job j to machine i. 

Consider a machine i such that < i < £,t+6- The set of jobs J'i which the sorting procedure 
allocated to machine i, was scheduled on a machine i' in S' (possibly swapping pairs of jobs in common 
mini-classes). Let j be a job such that S*{j) = i. In order to prove that pj > 2^'"'' it suffices to show that 
i' > Assume by contradiction that i' < S,t+4- We have Wf < W^^^^ < 2'^*+"' which holds since the 

S 

works in S are monotonically non-decreasing and Wf^_^^ < W^^^^ = 2 = 2^'+*. On the other hand, 
wf > Wg^g > 2^'+^-^ By (© and LemmafUl we have Wf* < Wf • (1 + e) < Wf{l + 12e)(l + e) < 



12 



^^^(l + Ue). Therefore, using (O we get 



< 777 by 



< 



1 - Ue 
2W^-?(1 + Ue) 



1 - Ue 
1 1 + Ue 



1 - Ue 



32 



contradicting r > 5. 



4 A dynamic programming for computing the best highly structured solu- 
tion 

In this section we show how to compute the optimal structured schedule S*. Our algorithm will use a 
dynamic programming procedure which is based on a shortest path (or an optimal bottleneck path) in a 
directed layered graph G = {V, E) with weights on its vertices. 

We will define a layered graph, in which the algorithm computes a path corresponding to an optimal 
solution with respect to a given goal function. Each layer of an index 1,2, ... ,m corresponds to a machine, 
and each vertex in one of these layers encodes a set of jobs which were scheduled prior to the current 
machine, and a set of jobs which were scheduled up to and including the current machine. The difference 
between these sets easily reveals the work of the current machine, and allows us to restrict the paths in the 
graph to schedules in which the works are monotonically non-decreasing. Given the work Wi of the current 
machine i, the weight of the vertex is the load of this machine Li, or f{Li) for a well-behaved function /. 
The edges between layers correspond to compatibility conditions which in particular enforce the condition 
that the works of machines are monotonically non-decreasing. The order of the layers is according to the 
speeds of the machines, that is, machines with higher speeds have a higher index of their layers, and subsets 
of machines with a common speed are ordered according to a fixed ordering of the machines. The graph 
which we will use allows us to find any structured schedule and maybe additional schedules. The schedule 
which will be found will be at least as good as the structured schedule whose existence we proved in the 
previous section. To distinguish between several optimal solutions, and to prioritize the possible outputs, 
we number all vertices of each layer with distinct integers, and we always search for paths whose reverse 
sequence of numbers along the path (that is, the sequence of vertices given from the end of the path towards 
the beginning of the path) is minimal (lexicographically) out of paths which have an optimal cost with 
respect to our goal function. This property allows us to assume that there exists a total order over the paths 
in the graph, and the algorithm always outputs the minimal path (according to this order) which is optimal 
in the current scenario. 

The graph G will encode in each layer all possible short histories of the magnitude vectors (which we 
call short magnitude vectors). There will be a starting vertex s, also seen as the layer of vertices of index 
0, and an end vertex t, also seen as the layer of vertices of index m + 1, and we always look for a path 
in G from s to t. Thus, V consists of m regular layers denoted as 1, 2, . . . , m (one for each machine) and 
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two additional layers and m + I. For every possible structured schedule, there will be an s — t path 
corresponding to it (and possibly additional s — t paths corresponding to other feasible schedules). 

A short magnitude vector ip = {ip^, ip^jip"^, . . . , ip^) for machine i is a vector consisting of seven con- 
secutive distinct values in a magnitude vector a (that is, there exists 1 < ^ < r(a) such that = vt+ri-b 
for 77 = 0, 1, 2, ... , 6). If this vector is associated with machine i, then = If = vt+b for some 
value of t, then if) = (f^, J/f+i, . . . , vt+&)- If the magnitude ^Jj^ is the largest magnitude in a, we will let 0^ 
be the fictitious value 00. Similarly, if ifj^ is one of the smallest five values in a, we add —00 as the first 
components of ij). We say that a short magnitude vector ^ is quasi-consistent with a schedule 5 is consists 
of six consecutive distinct values of a magnitude vector a such that S is quasi-consistent with a. 

Other than the entries which are —00 or cx), a short magnitude vector must be such that it can be a part 
of a magnitude vector. Thus, > 1/;^ -|- r -I- 1 for 77 = 0, . . . , 5. In addition, we define a list of allowed 
finite components. 

Lemma 16 For j E J let pj = [log2 Pj]- Then, for every possible subset J' C J of jobs whose total size is 
W, we have 

j€J k=0 

Proof. Let j G J' be a maximum indexed job in J'. Then, W > pj and W < j ■ pj < n ■ pj. Therefore, 
Pj < riog2 VF] < \log2{n ■ Pj)] = [log2 n + log2Pj] < [log2ra] + Pj and the claim holds. ■ 

Corollary 17 The number of possibilities of short magnitude vectors is 0{n^). 

Proof. The set of different values for each component in a short magnitude vector is 

|"log2 n]+l 

U U fe + A:}U{-oo,oo} 

jGJ fc=-l 

since we are only interested in magnitude vectors which are quasi-consistent with some schedule. Therefore, 
there are at most (n • (log2 n + 4) + 2)^ = 0{n^) different short magnitude vectors. ■ 

Next, we define the set A(^) of active mega-classes for a short magnitude vector if). A mega-class 
k belongs to A(V') if there exists a value of r/ = 0, 1, . . . , 6 such that \k — < r, and an alternative 
mega-class 7/;''+^ — r — 1 (and perhaps a smaller alternative mega-class consisting of a single alternative job) 
belongs to A{ip) if it is an alternative mega-class consisting of alternative jobs of mega-classes + r + 
1, . . . , 1^"^+^ - r - 1 for values of 77 = 0, 1, . . . , 5 for which V""^^ -i/j" >2r + 2. 

The motivation for this definition of ^(V") is that if machine i has a short magnitude vector ip, then all 
jobs of size at most 2^°'"^ are scheduled on machines with magnitude at most in any structured schedule 
that is quasi-consistent with ijj, i.e., strictly before machine i (since Oj = tj/^). Moreover, all jobs of size 
more than 2'^6+r+i 

are scheduled on machines with magnitude at least ip^ in any structured schedule that is 
quasi-consistent with ijj, i.e., after machine i. Thus the only relevant (alternative) mega-classes for machine 
i are the ones described above. These properties will be enforced by the structure of the graph. Moreover, 
given a set of consecutive mega-classes it can be decided to convert the jobs of these mega-classes into 
alternative jobs, and this can only happen if no jobs of these mega-classes were already scheduled. Once it 
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is decided, this decision is irrevocable and future sets of consecutive mega-classes which are converted into 
alternative jobs will be disjoint. 

A status vector of a short magnitude vector ijj consists of a component for each mini-class which belongs 
to a mega-class in A{'il)). This component represents the number of jobs (or alternative jobs if this is an 
alternative mini-class) which were already scheduled (recall that in a structured schedule we always schedule 
these jobs sorted by their sizes (with a fixed tie-breaking policy), and therefore the number of jobs which 
were scheduled uniquely identifies which jobs these are). 

Lemma 18 The number of status vectors for one specific short magnitude vector ip is ©(n^^^^'"'*'^)'^^^)^). 
Therefore, overall there are 0{n^^^^'^'^^^^^'^^^) status vectors. 

Proof. The claim holds since every component in the status vector is an integer in [0,n], the number of 
mini-classes in a mega-class is [logi+g 2] = A, and there are at most 12 alternative mega-classes in A{ip). 
m 

Consider a pair of ordered pairs {ip, u) and {ip', u') where u and u' are status vectors of the short mag- 
nitude vectors ip and V'', respectively. We say that such a pair is compatible if one of the following cases 
hold. 

1. If ^ = V^' and every component in u is at most its corresponding component in u' . 

2. If for all r/ = 1, 2, . . . , 6, V''' = (V'')^'^' and every component in n corresponding to a mini-class {k,t) 
(such that mega-class k is in A{il)')) is at most its corresponding component in u' . Moreover, every 
component in u' which corresponds to a mini-class (A;, t) such that k ^ A^if)) is zero. Informally, jobs 
of such zero components in u' are too large for ip. 

If ('0, u) and (-0', u') are compatible, then their difference defines a set of jobs which can be scheduled 
on a machine. This set of jobs J((V', ^0) defined as follows. The set J{{ip, u), {tJ/ ,u')) will 
contain all remaining jobs of mini-classes which have corresponding components in u but not in u' (these 
are the last jobs of each mini-class which are not scheduled yet, according to the information encoded in u). 
Informally, such jobs are too small for ip' and must be assigned immediately. For every mini-class which has 
components in both u and u', the number of jobs of this mini-class in J{{tp, u), {ip',u')) is the difference 
between these components (these are the next jobs in each mini-class). We denote by W{{iIj,u), {ip' ,u')) 
the total size of jobs in J{{'ip, u), {t/j' ,u')). 

The set of vertices of layer i (for i = 1, 2, . . . , m) is the set of compatible pairs {ip, u) and (t/;', u'). 
Thus such a vertex corresponds to ((V', u), {ijj',u')). The meaning of such a pair is to assign the jobs of their 
difference to machine i (and thus the work of i would be exactly W{{ip, u), {ip' , u'))), where -0 is the short 
magnitude vector of machine i, and ijj' is the short magnitude vector of machine i + 1. 

The weight of such a vertex in layer i is defined as '^((^'"^'(^ jf we are solving the minimum 
makespan problem or the problem of maximizing the minimum load. If we are interested in the problem of 
minimizing fi^i) ^ well-behaved function /, then the weight of the vertex is ,u )) y 

The vertices s,t do not have weights. 

A vertex in layer i (for 1 < i < m — I) corresponding to {{ip, u), {ip' , u')) is adjacent to a vertex in 
layer i + 1 corresponding to {{ip', u'), {ip" , u")) if and only if W{((ip, u), {ip' , u')) < W{{ip', u'), {ip" ,u")). 
There are no other edges between these layers, that is, there can be no edge from {{ip^,u^), {ip'^,u'^)) to 
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{ip'^, n'^)) in consecutive layers if (V''^, n'^) / {ip'^,u'^). The vertex s of layer is adjacent to all 
vertices in layer 1 corresponding to {{tp, u), {ip' ,u')) such that all components of the status vector u are zero, 
and = —oo. The vertices of layer m which are adjacent to t (of layer m+ 1) are the ones corresponding to 
((V', u), {ip',u')) such that ip'^ = oo, and for every mini-class whose mega-class is in A{ip') the component 
in u' is exactly the number of jobs in this mini-class (also for an alternative mini-class). The topology of the 
graph G depends only on the set of jobs and their sizes, and on the number of machines (and not on their 
speeds). Only the weights depend on the exact problem and on the speeds of the machines. 

We observe that an s—t path in the graph gives immediately a schedule, since each vertex {{ijj,u),{ijj' ,u')) 
in the graph defines a specific set of jobs allocated to the machine with index equal to the index of its layer, 
whose total size is exactly W{{'i/j, u), {tp', u')). Moreover, every s — t path defines a partition of the job set, 
and every such solution satisfies that the works of the machines are monotonically non-decreasing in the 
index of the machine. We also observe that every structured solution corresponds to (at least) one s — t path 
in the graph G. 

Using this graph, we compute a label for each vertex. This label is equal to the cost (or value) of the 
partial solution defined by the best path from s to this vertex. Moreover, we compute a pointer tt to the 
previous vertex on this best path from s. If there are several possibilities for best paths (ending at the same 
vertex) vr is defined to be the minimum index of the vertex satisfying these conditions according to the 
numbering of vertices in each layer. We next define the notion of a best path for each of the objectives 
considered in this paper. For the problem of minimizing the makespan, a best path is one that minimizes the 
maximum weight of a vertex along the path. For the problem of maximizing the minimum load, a best path 
is one that maximizes the minimum weight of a vertex along the path. Finally for the problem of minimizing 
Si^i fi^i) where / is a well-behaved function, a best path is a path of minimum total weight of its vertices. 



5 Monotonicity proof 

Our monotonicity proofs are based on analysis of a scenario where machine 7 changes its speed. We will 
assume that every machine 7' 7^ 7 has a fixed speed of s^i while machine 7 has two possible speeds s^ 
and s'^. We sometimes consider additional speeds between s-y and s'^. In the next two lemmas si, . . . , 
denotes a sorted hst of machines speeds. 

Lemma 19 Consider two executions of the algorithm, both with respect to minimizing Y^^i fi^i) where 
f is a well-behaved function (with a common function f), where the sorted order of machines is 1,2,..., m, 
each with its own set of speeds, resulting in the two schedules Si and S2 found by the paths Pi and P2. The 
two sets of speeds are defined as follows. For every i' ^ i the speed of i' is Si' in both sets, and the speed 
ofi is ai and 02, respectively, such that Sj-i < ci < 02 < Sj+i (where sq = and Sm+i = 00). Then, 

Proof. For a schedule S denote by COST5,cost|^ the costs of schedule S using the speeds ai and 02 
for machine i, respectively. Recall that the graph G remains the same in the two executions. Since the 
path Pi could have been found by the algorithm when it computes P2 and vice versa, COSt|^^ > COST^^^, 
and COST52 > COST^j, which gives COST'^^ - COSTs^ > COST'^^ — COSTg^. Assume by contradiction 
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Since cj2 > ai and wf^ < Wf\ we find WfH^ - — ) > W^H^ - Rearranging tlie last 

14/^1 w^'^ W^'^ w^^ 
inequality gives —^^ — I — ^ > Since i does not change its position in the sorted order of 

machines, we have COSt'^^ - COSTsi = /(^) - /(^) and COST^^ - cosTs^ = /(^) - /(^), 

and so we find using COSt'^^ - COSTs, > COSt'^^ - COSTs^ that /(■^) - f{^) > /(^) - /(^)- 

Using ai < a2 and W^"" < Wp we have ^ < < By convexity, we find / (-^)+/(-^) > 

/(^ ) + /(^ + ^ - ^^)- Using'strict mlotomaty of /, /(^ + ^ - #) > /(%), 
which is a contradiction. ■ 

The next lemma is used in the case that the speed of a machine changes. We will split the process of 
changing the speed into steps, and one type of step will be swapping the positions with another machine of 
the same speed. Therefore, we note the following. 



Lemma 20 Consider two executions of the algorithm, both with respect to the same objective function, each 
with the same set of speeds si, S2, ■ ■ ■ , Sm where Si = Sj+i, where the sorted order of machines is given by 
increasing indices in the first execution and the order obtained by swapping the positions of machines i + 1 
in the second execution, resulting in the two schedules Si and 82- Denote by wi the work of machine i in 

S S 

the schedule Si (that is, ooi = and by ijJ2 the work of the same machine in S2 (that is, uj2 = l^j+ij- 

Then, coi < 002- 



Proof. Since the set of optimal solutions for the two inputs is exactly the same, so is the set of optimal 
paths in the graph. Since our algorithm always outputs lexicographic minimal optimal path, we conclude 
that = ^2. The claim holds because the solutions obtained as paths in the graph have monotonically 
non-decreasing works of machines. ■ 

Theorem 21 The approximation scheme for minimizing X^i^i fi-^i) where f is a well-behaved function is 
a monotone PTAS. The approximation scheme for minimizing the £p norm of the vector of machine loads 
(obtained by running the algorithm with f{x) = x^) is a monotone PTAS even if p is a part of the input. 

Proof. Let S be an optimal solution, then by Theorem[T5l there is a structured solution S* such that for every 
i we have W^* < (1 + 14e) • , and thus Lf * < Lf • (1 + 14e), and therefore the cost of S* as a solution to 
our problem is at most E^i/l^f) < TT=i f i^l ■ {I + Ue)) < (l + 0(l)e) ^™ ^ /(Lf) where the first 
inequality holds by monotonicity of /, and the second inequality by the property of / that ifx < (l + e)y 
then f{x) < (1 + 0{l)e)f{y). The schedule given by the algorithm as output has a cost which is no larger 
than the cost of S*. Note that the approximation ratio of S* for the problem of minimizing the Ip norm of 

the vector of machine loads is at most 1 + 14e since {YT=i i^fYf^^ < (1 + • (E^i {LiY)^^^- 

To prove the monotonicity, consider a machine i which increases its speed from Sj to s[. We split the 
process of increasing the speed of a given machine into two types of events. The first type are time intervals 
in which the position of this machine in the sorted order of the machines does not change. The second 
type are points in time when the speed is fixed, but the machine swaps its location with the next machine 
in the list of machines sorted by speed. There can be multiple such time intervals and points in time, and 
it is sufficient to consider one event of each type, thus we consider two cases. The case where machine i 
increases its speed, < Sj+i, and machine i does not change its position in the sorted list of machines, and 
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the case Sj = Sj+i, where the only change is that these two machines swap their relative order. For the first 
case, the claim follows by Lemma [19] For the second case, the claim follows by Lemma |20l ■ 

Theorem 22 The approximation scheme for maximizing miiijgA/ L-i is a monotone PTAS. 

Proof. Let 5* be an optimal solution, then by Theorem [151 there is a structured solution 5** such that for 
every i we have W^* > (1 — 14e) • , and thus Lf * > if • (1 — 14e), and therefore the value of S* is 
at least 1 — 14e times the value of 5. The schedule given by the algorithm as output has a value which is no 
smaller than the value of S* . 

To prove the monotonicity, consider a machine i which increases its speed from Si to s'^. Consider the 
solution Si obtained by the algorithm for the case where the speed of i is Sj. Let Ci be the value of Si 
(computed for the set of speeds where the speed of i is Si). We split the process of increasing the speed of 
a given machine into two periods where the first period is split further into two types of events. In the first 
period, the speed of i is at most a, where a is the maximum speed for which the value of the solution Si is 
exactly Ci (possibly swapping the contents of machines if machine i changes its position in the sorted list of 
machines according to the sorting done by the algorithm). Note that a is well-defined, that is, the maximum 
exists. If fj = Sj, we say that this period is empty. If cr > s^, we set a = s[. Therefore, during the first 
period the speed of i is in (sj, a]. If u = s[, then the second period is empty, and otherwise the speed of i is 
in (cT, s ■] in this period. For the first period, the first type of events are time intervals in which the position of 
this machine in the sorted order of the machines does not change. The second type are points in time when 
the speed is fixed, but the machine swaps its location with the next machine in the list of machines sorted by 
speed. 

We prove that for every speed in [sjjfj], the solution is returned by the algorithm. First, we show 
that the value of an optimal solution remains Ci. The value of an optimal solution cannot increase when i 
increases its speed, so by definition Si remains an optimal solution. Moreover, when i increases its speed 
in the first period, the set of optimal solutions is a subset of the set of optimal solutions when the speed of 
i is Si (even if locations of machines are swapped). Therefore, the algorithm outputs for every speed in 
the first period. Thus, for time intervals in which the position of i in the sorted list of machines is fixed, the 
work of i is exactly the same, and in events in which machine i swaps its position with another machine, the 
work of i cannot decrease by Lemma |20l In the case o" = s ■ we are done. Otherwise, we assume that there 
are no further machines of speed a which appear later than i in the ordering of the machines (possibly by 
adding events of the second type for the first period). 

Next consider the case where a < s'.-. Denote by W the work of i in the solution Si where the speed 
of i is a. Recall that for this speed of i, the value of the optimal solution (i.e., of Si) is exactly Ci. We 
prove that ^ = Ci. Assume by contradiction that the claim does not hold (that is we assume that ^ > Ci, 
as otherwise the value of Si in this case is strictly smaller than Ci contradicting the definition of cr). Let 
(Ti > (7 be such that ai < ^ and ui is smaller than the speed of the next machine after i in the sorted list 
of machines, if such a machine exists. Then, the value of Si for the speed ai of i remains Ci contradicting 
the maximality of a. Let C2 be the value of an optimal solution 52 found by the algorithm where the speed 
of i is Then, C2 > Ci • since otherwise Si is a strictly better solution for speed s[ of i, because even 

s' 

if machines swap locations the machine in every position is faster by no more than . Denote by W the 
work of iin S2- We have W > C2 ■ s'^ > Ci ■ a = W, and the claim follows. ■ 

The proof of the next theorem is similar to the proof of Theorem |22l and it is given in Appendix IA.2I 
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Theorem 23 The approximation scheme for minimizing maxjgjvf Li is a monotone PTAS. 

6 Computing the payments 

Archer and Tardos [5] defined a payment scheme which can be applied for any monotone scheduling algo- 
rithm to create a truthful mechanism. Denote the payment to agent i by Pj. We briefly repeat the definition 
of Pi. Let b-i denote the vector of bids, not including agent i. We write b (the complete bid vector) also as 
{b_i, bi). Then the payment function for agent i is defined as 



where Wi{b-.i, bi) is the work (total size of jobs) allocated to machine i given the bid vector b and the hi are 
arbitrary functions (Theorem 4.2 in |5]). 

In order to compute the payments, we need to calculate the integral in ([5]). Recall that the bid of an agent 
represents its claimed cost for processing one unit of work, which can be seen as the inverse of the speed 
of its machine. For a given set of bids . . . , bm), calculating the integral for agent i requires us to know 
what its work would be for every possible bid b of this agent, i.e., for the bids b) for b € (0, oo). First, 
we partition the possible bids into intervals in which the position in the ordered set of machines (that is, its 
layer in the graph) of machine i remains constant. Consider the set {0, oo} U and denote its 

elements by = ci < • • • < Cm' = oo (m' < m + 1), then the intervals to consider are {cj, Cj+i) for 



For each vertex v in layer i, we compute a function -F^(6) which is the objective function value of the 
best path which traverses this vertex, as a function of the bid of machine i. Recall that the algorithm outputs 
the minimum or maximum (over all vertices of the layer) of these functions depending on the objective 
function. 

Claim 24 For each vertex v and every bid interval {cj, Cj+i), F^(b) is a piecewise linear continuous func- 
tion with a polynomial number of pieces. 

Proof. In layer i, the weight of vertex v which represents the compatible pair ((■(/', u), {ip' ,u')) is the constant 
W{{tlj,u), ('0', u')) divided by Si, where Si = l/bi. Note that the pair represented by v also specifies the set 
of jobs assigned to machines before machine i, and the set assigned after i. Due to the tie breaking done in 
the dynamic program, and the fact that only the speed of machine i changes, this means that the identity of 
the best path which passes through v does not depend on 6 (only its objective value does). 

For the makespan and the maximizing the minimum load problems, the objective value of a path is the 
maximum (minimum, respectively) weight of a vertex along the path. Hence, as bi increases from Cj to 
Cj+i, the only change that can happen is that the weight of vertex v starts having the maximum weight 
along the fixed best path (for the makespan objective) or stops having the minimum weight (for the covering 
objective). Therefore, Fy{b) has at most two pieces, where for one piece machine i is a bottleneck machine 
(that is, a machine whose load equals the objective function value of the solution) and for the other it is not. 
If i is the bottleneck, Fy{b) = b • W{{'ip, u), {ip' , u')), else Fy{b) is constant. 

For the minimization of Yl^i fi^i) for a well-behaved function /, the objective value of a path is the 
total weight of its vertices. Here, Fy{b) is a constant plus /(6 • W{(ip, u), ,u')) (where the constant is the 




(5) 



j = l,...,m' -1. 
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total weight of the other vertices along the best path which traverses v). Hence by using the approximated 
piecewise-linear convex monotonically increasing function of / instead of / itself the claim follows since it 
is sufficient to consider such an approximated function with pieces ending at integer powers of (1 + e) (and 
thus with polynomially many such pieces). ■ 

Claim |24] implies that the number of intersection points between any pair of functions Fu{h)) is 

also polynomial. Thus we can compute all of these points in polynomial time, and determine which points 
lie inside the interval (cj, Cj+i). Moreover, we can also determine which schedule our mechanism uses for 
each intersection point by running the PTAS for each point, including Cj (if Cj > 0) and Cj+i (if Cj+i < oo). 
After removing duplicates, this gives us a list of intersection points with associated schedules and works. 

Remark 25 The replacement of f with the convex monotonically increasing piecewise-linear approxima- 
tion of f is crucial. Without it, computing the value ofbi in which one solution becomes better than another 
solution involves computation of an exact solution of equations involving convex functions ( this cannot be 
done even for the case where f{x) = x^). However, for piecewise-linear functions this can be done effi- 
ciently. 

It is now straightforward to determine the schedule used for any possible bid b, and from that the work 
for any bid, as follows. Note that the schedule chosen does not change between any pair of consecutive 
intersection points by construction. Thus the work remains constant between any such pair. If the schedule 
used is the same at both endpoints, the work in between is given by this schedule. If two different schedules 
are used, then in the entire open interval between the pair, the schedule is used which gives the best value 
for the objective function. This can be determined by running the PTAS for one point inside this interval. 
Thus we can find the exact value of the integral in ^ (without rounding the speeds of the machines). 

A Omitted proofs 

A. 1 Proof of the second part of Claim [1] 

Consider a schedule S with makespan M and cover C. Call a pair of machines i, j reversed HI < i < j < 
m and Wf > W'^. We show that removing a consecutive reversed pair (that is, j = i + 1) by swapping 
the sets of jobs assigned to them from any schedule S does not increase the makespan or decrease the cover, 
which implies the claim (since after a finite number of such steps there will no longer be reversed pairs). 
Let S' be the schedule resulting from swapping the two job sets of machines In S', machine j gets more 
work, but the load remains at most M: we have /sj = Wf /sj < Wf /si < M. Machine i gets less 
work, but the cover remains at least C : we have Wf / Si = /si> / Sj > C. 

A.2 Proof of Theorem |23] 

Let S be an optimal solution, then by Theorem [15] there is a structured solution S* such that for every i we 
have wf* < (1 + 14e) • Wf, and thus Lf * < Lf • (1 + 14e), and therefore the makespan of S* is at most 
1 + 14e times the makespan of S. The schedule given by the algorithm as output has a makespan which is 
no larger than the makespan of S*. 

To prove the monotonicity, consider a machine i which decreases its speed from Si to sJ. Consider the 
solution Si obtained by the algorithm for the case where the speed of i is Sj. Let Ci be the makespan of Si 
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(computed for the set of speeds where the speed of i is Si). We spht the process of decreasing the speed of 
a given machine into two periods where the first period is split further into two types of events. In the first 
period, the speed of i is at least a, where a is the minimum speed for which the makespan of the solution 
Si is exactly Ci (possibly swapping the contents of machines if machine i changes its position in the sorted 
list of machines according to the sorting done by the algorithm). Note that a is well-defined, that is, the 
minimum exists. If a = Si, we say that this period is empty. If cr < s-, we set a = s-. Therefore, during the 
first period the speed of i is in [a, Sj). If cr = s[, the second period is empty, otherwise the speed of i is in 
[s^, cr). For the first period, the first type of events are time intervals in which the position of this machine 
in the sorted order of the machines does not change. The second type are points in time when the speed is 
fixed, but the machine swaps its location with the previous machine in the list of machines sorted by speed. 

We prove that for every speed in [cr, Sj], the solution 5*1 is returned by the algorithm. First, we show 
that the makespan of an optimal solution remains Ci. The makespan of an optimal solution cannot decrease 
when i decreases its speed, and by definition Si remains an optimal solution. Moreover, when i decreases 
its speed in the first period, the set of optimal solutions is a subset of the set of optimal solutions when the 
speed of i is Si (even if locations of machines are swapped). Therefore, the algorithm outputs Si for every 
speed in the first period. Thus, for time intervals in which the position of i in the sorted list of machines 
is fixed, the work of i is exactly the same, and in events in which machine i swaps its position with the 
previous machine, the work of i cannot increase by Lemma [20l In the case cr = s • we are done. Otherwise, 
we assume that there are no further machines of speed cr which appear earlier than i in the ordering of the 
machines (possibly by adding events of the second type for the first period). 

Next consider the case where a > s'^. Denote by W the work of i in the solution 5*1 where the speed of 
i is cr. Recall that for this speed of i, the makespan of the optimal solution (i.e., of ^i) is exactly Ci. We 
prove that — = Ci. Assume by contradiction that the claim does not hold (that is we assume that ^ <Ci, 
as otherwise the makespan of in this case is strictly larger than Ci contradicting the definition of cr). 
Let cri < cr be such that cri > ^ and cri is larger than the speed of the previous machine before i in the 
sorted list of machines, if such a machine exists. Then, the makespan of 5i for the speed cri of i remains Ci 
contradicting the minimality of cr. Let C2 be the makespan of an optimal solution S2 found by the algorithm 
where the speed of i is s[. Then, C2 < Ci • since otherwise 5*1 is a strictly better solution for speed 
of i, because even if machines swap locations the machine in every position is slower by no more than 
Denote by W the work of i in 82- We have W' < C2 ■ s[ < Ci ■ a = W, and the claim follows. 
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